Molecular gyroscopes and biological effects of weak ELF magnetic fields 



V.N. Binhi* 

General Physics Institute Russian Academy of Sciences, 38 Vavilova St., 119991 Moscow, Russia 

A.V. Savin* 

Institute for Physics and Technology, 13/7 Prechistenka St., 119034 Moscow, Russia 

Extremely-low-frequency magnetic fields are known to affect biological systems. In many cases, biological 
effects display 'windows' in biologically effective parameters of the magnetic fields: most dramatic is the fact 
that relatively intense magnetic fields sometimes do not cause appreciable effect, while smaller fields of the 
order of 10-100 /iT do. Linear resonant physical processes do not explain frequency windows in this case. 
Amplitude window phenomena suggest a nonlinear physical mechanism. Such a nonlinear mechanism has 
been proposed recently to explain those 'windows' . It considers quantum-interference effects on protein-bound 
substrate ions. Magnetic fields cause an interference of ion quantum states and change the probability of ion- 
protein dissociation. This ion-interference mechanism predicts specific magnetic-field frequency and amplitude 
windows within which biological effects occur. It agrees with a lot of experiments. However, according to 
the mechanism, the lifetime of ion quantum states within a protein cavity should be of unrealistic value, 
more than 0.01 s for frequency band 10-100 Hz. In this paper, a biophysical mechanism has been proposed that 
(i) retains the attractive features of the ion interference mechanism, i.e., predicts physical characteristics that 
might be experimentally examined and (ii) uses the principles of gyroscopic motion and removes the necessity 
to postulate large lifetimes. The mechanism considers dynamics of the density matrix of the molecular groups, 
which are attached to the walls of protein cavities by two covalent bonds, i.e., molecular gyroscopes. Numerical 
computations have shown almost free rotations of the molecular gyros. The relaxation time due to van der Waals 
forces was about 0.01 s for the cavity size of 28 angstroms. 



I. INTRODUCTION 

Weak static and extremely-low-frequency (ELF) magnetic 
fields (MFs) can affect living things: cells, tissues, physiolog- 
ical systems, and whole organisms [1, 2, 3]. In many cases bi- 
ological effects of weak MF feature resonance-like multipeak 
behavior. Multipeak responses or magnetobiological spectra 
may appear with varying the frequency or amplitude of AC 
MF [4] and the magnitude of DC MF [5]. Usually, the term 
'windows' is used for the peaks of the spectra. 

Amplitude 'windows', see Fig. 1, specify nonlinearity of 
the transduction mechanisms involved in magnetobiological 
effects. This is confirmed more by the fact that magnetic 
noise simultaneously superimposed on a regular magnetic sig- 
nal suppresses biological effect of that signal [6, 7, 8, 9]. 

A nonlinear mechanism based on quantum interference 
has been developed in [10] to explain unusual ELF MF fre- 
quency and amplitude dependencies of magnetobiological ef- 
fects (MBEs). The mechanism elaborates the interference of 
ions bound within proteins. According to this mechanism, su- 
perposition of the ion states forms a non-uniform pattern of 
the probability density of ion. This pattern consists of a row 
of more or less dense segments occurring due to the inter- 
ference between quantum states of ions in a protein binding 
cavity. In a DC MF the pattern rotates with the cyclotron fre- 
quency. Exposure to a time-varying MF of specific parameters 
retards the rotation of the pattern and facilitates escape of the 
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ion from the cavity. This escape might influence equilibrium 
of biochemical reactions to ultimately result in a biological 
effect. 

Biologically effective parameters of AC-DC magnetic 
fields depend on the charge-to-mass ratio of the ion in ques- 
tion. The closed formula is derived for 'magnetic' part P of 
ion-protein dissociation probability. Predictions based on this 
formula reveal good agreement with experimental results in- 
volving calcium, magnesium, potassium, hydrogen and other 
ions of as molecular targets for MF. The theory describes mul- 
tipeak frequency and amplitude spectra of MBEs involving 
ions of Ca 2+ , Mg 2+ , and H + as molecular targets for AC-DC 
MFs [10]. 

The interference mechanism is surprisingly effective in ret- 
rospectively predicting results of existing experiments con- 
ducted under the following defined MF conditions: parallel 
AC-DC and pulsed MFs [10, 11], 'null' and static MFs [13], 
and various MFs with a slow rotation of a biological system 
[12]. As an example, Fig. 1 demonstrates the comparison of 
experimental data, in parallel AC-DC MFs, on MBEs involv- 
ing fixed and rotating proteins, and calculated curve (dash 
line). 

The good consistency between theoretical calculations and 
many experiments indicates that what underlies magnetobio- 
logical effects is most likely an interference phenomenon. 

According to the interference mechanism, the relation 
should be valid r _1 fi c > 1, where is the lifetime of 
ion quantum states within a bound cavity and Sl c is a cy- 
clotron frequency of an ion in the geomagnetic field, usually 
10-100 Hz. The postulate therefore has to be made that ion 
quantum states, more exactly their angular modes, live more 
than 0.01 s within the cavity. However it is in contradic- 
tion with our common knowledge that such states might live 
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FIG. 1: Experimental evidence [17, 18, 19, 20, 21, 22, 23, 24, 25, 26] for MBEs in a uniaxial MF ([18] noted a weak perpendicular component 
of a DC MF). Theoretical amplitude spectra: a dash line was derived for fixed ion-protein complexes (factor a = 1 was not shown) and also 
for rotating ion-protein complexes, see details in [12]. Solid line represents the function (18) derived for gyro interference. 



only 10~ 12 -10~ 10 s because of the thermalizing interaction of 
ion with cavity walls. On the other hand, the weak AC MF, 
Kl c -c fc B T, is commonly believed to be unable to contribute 
into thermally driven (bio)chemical reactions (so-called kT- 
problem). 

To overcome the problem, we note that there is a specific 
mechanism that provides relatively large lifetime of the an- 
gular modes. Consider a dipole molecular group that are at- 
tached within the cavity to its walls in two points, i.e. by two 
covalent bonds, thus forming a group that may rotate inside 
the cavity without contact with walls. Such a construction is 
referred to as gyroscope. In the case, it is a molecular gy- 
roscope. Of importance is the fact that thermal oscillations 
of that covalent bonds, or gyroscope's supports, make only 
zero torque about the axis of rotation. This leads to rela- 
tively slow thermalization of a gyroscopic degree of freedom. 
Relaxation is mainly due to van der Waals interaction with 
thermalizing walls. As far as the interaction potential, the 
Lennard-Jones potential, decreases as r~ 6 and walls' inner 
surface grows as r 2 , the overall van der Waals contribution 
varies approximately as r~ 4 . That is, relaxation quickly di- 
minishes with the cavity size to grow. Computations show 



almost free rotations (thermalization time 0.01 s) of a molecu- 
lar gyro within the cavity of 28 angstroms size. This is enough 
for the ion interference mechanism to display itself. Probably, 
such roomy cavities are formed by ensembles of a few protein 
globules, between them, or within some enzymes that unfold 
DNA double-helix. 



n. MOLECULAR GYROSCOPE 

A long lifetime of angular modes is the sole serious ideal- 
ization underlying the mechanism of ion interference. This 
idealization would be hard to substantiate with the ion-in- 
protein-capsule model. One would have to assume that the 
ion forms bound states of the polaron type with capsule walls. 
In turn, justification of a large lifetime of polaron angular 
modes would require new idealizations. A 'vicious circle' 
occurs which one could not leave without having to substan- 
tially change the model itself. Thus, despite the obvious ad- 
vantages of the ion-in-capsule model, namely, simplicity and 
a high forecasting skill, we have to recognize its limitations 
and seek for other solutions. 
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FIG. 2: Forces, moments of forces, and angular momenta in rotation 
of a gyro. 



One of them hinges on the use of conservation laws in the 
dynamics of rotating solids. Rotation of a solid is described 
by the equation 



(1) 



where L is the angular momentum, K is the sum of torques 
acting on the solid. Consider for simplicity a symmetric gyro 
rotating around one of its main axes of inertia with a force F 
acting on its point of support, as shown in Fig. 2. The moment 
of this force about the shown axis is obviously zero. From 
equation (1) we have 

L = L + dL, dL = K dt = r x Fdt. 

Since K _L F, then dL _L F, i.e., the force caused an orthog- 
onal displacement of the axis of rotation. Also, the vector r is 
directed along the axis of rotation, therefore the vector dL is 
also orthogonal with Lo . 

Thus, a continuously acting force F causes a forced preces- 
sion of the gyro about the direction F with an angular velocity 
defined by the angle through which the gyro axis of rotation 
deviates per unit time, viz., 



n 



precession 



riL/Lrj 
dt 



K _ rF 



The length of vector r is defined by the gyro locking condi- 
tions. If point B is fixed, then the origin of r coincides with 
B. If point B is free, then the origin of r is on line AB and de- 
pends on the gyro parameters. For estimation, it is important 
that r has the order of magnitude of gyro length. 

Let the gyro be a model of a rigid molecule free to move 
and constrained by the thermal oscillations of one of the point 
of support (e.g. A) alone. We estimate the mean gyro axis 
deviation angle for a random force F causing chaotic oscilla- 
tions of its point of support. It should be noted that the gyro 
gravity energy ~ MgR is many orders of magnitude below its 
kinetic energy ~ L 2 /2I and the effects of gravity may be ne- 
glected. In the last formulas, M, R, and I are the gyro mass, 



size, and moment of inertia, and g is the acceleration due to 
gravity. 

The energy of natural gyro rotation is Eq = Lq/2I. The 
gyro energy including chaotic rotations is Eq + k B T. On the 
other hand, the mean energy with allowance for orthogonality 
of Lq and dL is 
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where brackets mean averaging over the ensemble. Then 
(d 2 L)/2I ~ k B T. Denoting the average deviation angle by 
a = y / (d 2 L}/L yields a 2 ~ 2Ik B T/L 2 . The smaller L 
the larger the random deviations of a molecule caused by ther- 
mal perturbations of its support. Such a support is the covalent 
bond with the body of protein molecule. Low bound estimates 
of Lq follow from the Heisenberg uncertainty principle which, 
for a complementary pair of noncommuting operators of an- 
gular variable cp and angular momentum C ~ d / dip, can be 
written as: 

ALAip ~ ft/2. 

Since Aip ~ it, then AL ~ h/2n; thus the angular momen- 
tum cannot be smaller than its uncertainty, i.e., Lo ~ h/2ir. 
Finally, we have 
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As can be seen deviations increase with the size of molecule; 
however, even for small molecules, the estimate of deviation 
is unrealistically large. It implies that, in lower rotation states, 
molecules will 'lay aside' in response to perturbation of their 
support and, consequently, the angular momentum will not 
be conserved. It should be noted that we are interested only in 
angular states with small quantum numbers. Otherwise the in- 
terference patterns to be discussed below become fine grained 
and are unlikely to be reflected in measured properties. 

Thus, in order to be immune to thermal displacements of 
supports, the gyro has to have its second support also fixed in 
the protein matrix. The configuration of a rotating solid with 
supports fixed in the rim is one of the types of a gyroscope, 
i.e., a device to measure angular displacements and velocities. 
What we consider is essentially a molecular gyro: a relatively 
large molecular group is placed in a protein cavity and its two 
edges form covalent bonds (supports) with the cavity walls. It 
is important to note that thermal oscillations of the supports 
produce only zero moments of forces about the natural group 
rotation axis. Therefore, the gyroscopic degree of freedom ip 
is not thermalized by the supports' oscillations. This does not 
imply that the energy of the gyro does not dissipate. Radiation 
damping or Lorentz friction force is neglected, because of its 
infinitesimal value. Below we examine at first the interference 
of the molecular gyro and then the damping due to wan der 
Waals forces. 
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FIG. 3: An amino glutaric acid molecule with potentially ionizing 
groups. The z axis is the main axis of inertia. Rotation of charges 
distributed over the molecule in a magnetic field leads to interference 
of its quantum angular states. 



We choose the vector potential in the form 

-±Hy, ±Hx, 



(5) 



With allowance for constraints (4), the velocity of a particle in 
spherical coordinates will be v = Rip, and the velocity vector 
in Cartesian coordinates is 



-Ripsm(tp), Ripcos(ip), 0) 



(6) 



Substituting this expression in equation (3), we obtain the La- 
grange function in spherical coordinates 



L = —— + —R^- q A . 



(7) 



Now, the generalized momentum is I = dL/dip, and the 
Hamilton function H = Ip — L is equal to 



III. INTERFERENCE OF THE MOLECULAR 
GYROSCOPE 

Rotations of large molecules is much slower a process than 
electron and oscillatory processes. Therefore, we think of the 
rotating molecular group as a rigid system of charged point 
masses — atoms and molecules with partially polarized chem- 
ical bonds. To illustrate, we point to molecules of amino 
acids which could be built into rather spacious protein cavi- 
ties forming chemical bonds at extreme ends of the molecule, 
thus forming a molecular gyroscope. Amino acids are links 
of polymeric protein macromolecules and also occur in a bio- 
plasm as free monomers. The general formula of amino acids 
is well known: 



R 



H 2 N+ - CH - COHCT , 

where R is a radical which differs one molecule from another. 
Polarities of the groups are shown in a water solution. By way 
of example, the radical of amino glutaric acid consists of three 
links — CH2 — CH2 — COOH, as shown in Fig. 3. Fixed on 
either side of a cavity, such a molecule, treated as a dynamic 
unit, has one degree of freedom — a polar angle ip, which 
simplifies analysis of its behavior in a magnetic field. 

For small velocities, the Lagrange function of one charge 
particle has the form 



Mv 2 



(3) 



where v is the particle velocity, and q is a charge. Let the 
magnetic field H = (0, 0, H) be directed along the z axis, 
and the particle be bounded by a holonomic constraint caus- 
ing its circumferential motion in the xy plane. In spherical 
coordinates, the constrains can be written in the form 



R = const. 



tt/2. 
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qA . 



(8) 



In the absence of electromagnetic field H = l 2 /2M R 2 , and it 
is obvious that I is the angular momentum of the particle. The 
Hamiltonian operator repeats (8) with the difference that here 
I is the angular momentum operator C = —ihd/dip. 

Let now a few particles rotate and, in a spherical system of 
coordinates, the constraints for particle i be 

ri = const. , 9i = const. 

Then, for a system of particles in a uniaxial magnetic field, the 
Lagrange function can be written following the derivation of 
formula (7) as 

j HQ -v 

L = 2^ + ~2c L P" ^2& A o(n,6i,<Pi) , (9) 



where 



/ = ]TAf 4 r 2 sin 2 (^), Q-^q^sin 2 ^) (10) 



is the moment of inertia, and 'charge moment of inertia' of 
the system about the axis of rotation. As can be seen, the 
Lagrange function of the system follows from the Lagrange 
function (7) after formal replacement of MR 2 with /, qR 2 
with Q, and qAo with the respective sum. Therefore, the 
Hamiltonian of the system immediately follows from equa- 
tion (8) after similar substitutions 



QHY 
2c / 



qiA (n,6i,(pi). 



(4) 



We assume further that the electric field is absent, i.e., let 
A Q = 0: 



2c J 
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In addition to C 2 /2I we find here two more operators. There 
are certain grounds to neglect the term proportional to squared 
H. From the ratios of coefficients at the terms quadratic and 
linear in H we obtain QH/Ach ~ 10~ 7 , where, for estima- 
tion purposes, we let Q ~ eR 2 , R ~ 10~ 7 cm, H ~ 1 G. 
Dropping this term we write the Hamiltonian in a convenient 
form 



(ii) 



The eigenfunctions and energies of the time-independent part 
of Hamiltonian (11) are 

l h 2 
— exp(imtp) 7 m = 0, ±1, e m = — m . 
2ir 21 



We now consider the ensemble of gyros that features a density 
operator a obeying the Liouville equation 



iha = Ha — ali. , a = 



(12) 



Some physical quantities, like the intensity of a spontaneous 
emission or the radiation reemitted by an ensemble, are known 
to linearly depend on the density matrix of the ensemble 



mm' 



The probability of biochemical reaction that we examine here 
is not a quantity of that sort. The reaction probability does 
not directly depend on the density matrix of the ensemble. 
It is rather the probability of the reaction of a gyro averaged 
over that ensemble. Therefore, at first we will find the density 
matrix <t^L of the oth gyro, then the reaction probability of 

(a) 

that gyro that non-linearly depends on <? m Li, and at last we 
will average the result over the gyro ensemble. 

Let the ensemble consist of gyros that appear with a con- 
stant rate at random moments of time. We assume the new 
gyros appear in a quantum state that is a superposition of the 
states close to the ground one, i.e., 



» 



,(0) = 



const , to, to ~ 1 
. to, to' ^ 1 



In the process of thermalization, the levels turn out to be popu- 
lated with the energies up to e m ~ k B T, i.e., with numbers up 
to to ~ i ^Ik B T ~ 10 3 for gyros with the inertia moments 
of the order of I ~ 10 35 g cm 2 . However, we are interested 
in the dynamics of the lowest states, that only could result in 
observable effects. 

In the representation of the eigenfunctions of 7i$ the den- 
sity matrix equation may be written from (11) and ( 12) as fol- 
lows 



"(rmm'^t"^mm.')^mm f ^ ^ ^ (VmjQ7m' ^ml^lm') 

(13) 



where 

w mm ' = ^( m2 - m ' 2 ) , V m i = -hu(t)mS ml . 

Phenomenological relaxation of the density matrix elements 
is taken into account, through the damping constants r mm < . 
Because of the relaxation the elements cr mm i of the lowest 
modes decrease while those of upper modes increase. As far 
as the stationery dynamics of a separate gyro is out of interest, 
we don't allow for the pumping upper modes, i.e., population 
redistribution into the states with large numbers m. 

Substitution of the above relations in (13) gives rise to the 
equation 

a = —Ya + ia [{in — m')u>{t) — u] , 

where indices to, to' are temporarily omitted for convenience. 
Along with notation 

g(t) = -T + if, f=(m- ro>(t) - u 

the equation takes the straightforward form & = g(t)cr. In the 
solution of that equation a = Cexp (J git) dtj, the constant 
C follows starting conditions. 

Let the MF possesses both DC and AC parts, then 

n\ 1 1 i w r>,n QHuc ,/ Hag 
Uj(t) = w g (l + h cos fit) , w g = — - — , h EE — . 

lie Hdc 
Now we separate constant and alternating parts in g(t): 

g(t) = —x + iztt cosilt , x ee r + iuj — i(m — to )w g , 

hi h' O 

zee (m-m')wg— = (to -to')— , fi' ee — . 

\ l M UJg 

The integral equals 

J git) dt = J (—x + izQ cos flt)dt = —xt + iz sin fit , 
hence 

a = cr(0)e^ 9{t) dt = a(0)e~ xt e tz sin nt 
= a(0)e~ Kt ^J„(z)e mnt . 

n 

Restoring indeces to, to', we arrive at the equation 

u mm' — u mm'\ VJ )^- 



Further, all the damping constants are assumed to equal T. 
With the notation 

(3 ee r + iui mm i — 7(777 — 777,')(jjg — infl , 

we rewrite the last equation in the form 
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that will be used later. 

Now we consider the probability density of a gyro to take 
an angular position tp, which is the only favorable position of 
the rotating group of the gyro to react with the active site on 
the wall 

p(t) = **(t, p)*(t, <p) = ± E c*m(t)e- lmiP E c ™' (*) e * m 

m m' 

- J_ \ ~" a -i(m-m')if 
mm' 

that is, 

p( t ) = ir E w(o)e- i(n, - m ' )v e-' 3t J Il ( 2mm ,). 

Z7T z — " 

mm'n 

i 



It is expedient to perform a sliding averaging in order to 
smooth out the relatively fast oscillations: they do not affect 
the active site that features character time constant r, i.e., 

i r t+T 

Pr(t) = jr- / p(t')dt' . 

Virtually, the factor exp(— fit) should be averaged: 




therefore 



Pr(t) = f E W(0)^^e-^— >e"^J„(w) • (14) 



Then, as in the ion interference model, we assume the reaction 
probability of a side group of the rotating molecule with the 
protein active site to be a non-linear function of the probability 
density (14). In the absence of whatever information on that 
function, it makes sense to consider quadratic non-linearity, 
since the linear term makes no contribution to that probability, 
see details in [10]. To find the reaction probability we will 
square (14) and take the average over the gyro ensemble. 

In the product p T (t)Pr{t) there are (i) complex conjugate 
terms, i.e., pairs with indeces n, m, m' and —n, m', m, which 
apparently do not oscillate, and (ii) fast-oscillating terms that 
we omit in view of the subsequent averaging. Omitting also 
immaterial numerical coefficient, we write 



Pr(*) 



,-2rt 



E kw(0)| 5 



sinh(/3r) 



3t 



J n i^mm' ) 



In this expression, the multiplier 



S= E kmm' (0)| 5 



sinh(/3r) 



0T 



contains the magnetic field dependence. 

Let a gyro appear in a moment of time t', then the reaction 
probability at time t equals 



u(t,t') 



Se -2T{t-t') ; t > t , 
, t<t' 



Assuming the moments of time t' to be distributed over the 
gyro ensemble in the interval (—0, 6) with a uniform density 
W (instead of a discrete distribution for wt^ in (12)), we find 
the mean probability P by proper integrating over the param- 
eter t'\ 

/e w g 
u(t, t') d£ = — — - . 
-e 21 



To link this value to an observable, e.g., a concentration of 
the reaction products, we write the kinetic equation for the 
number N of gyros per unit of tissue volume 



N 



w 



PN 



that gives TV = w/P = 2T/S in stationery conditions. Let So 
and TVo stand for corresponding quantities in the absence of 
an AC MF, i.e., at h' = 0. We would like to know the relative 
change p of the concentration of the reaction products under 
the AC MF influence. This is the relative number of gyros 
entering the reaction, i.e., 



Nq-N 
Nn 



1 



So 



(15) 



We now estimate values of S and p. The following notation 
will be used: 

/3t = ?7-M£, t] = Tt, £ = [u> mm < — (m — m')u) g — n£l]T . 
Then the expression for S takes the form 



E 

mm'n 



2 sinh rj + sin £ 2 



(m-m ) — 



(16) 

Since 77 is a constant, the frequency spectrum is defined 
mainly by the equation £ = 0, i.e., 

u m m' — uj g (m — to') — nfl = . 

For arbitrary small m, m! frequencies ui mm ' fall into the mi- 
crowave range. The effects of low-frequency MFs are defined 
by the interference of the levels to' = — m, when to mm i = 0. 
Then 

ujg (to — to ) + nil = , 
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from which we find 



0' = ^. 

n 



(17) 



The series over n in (16) converges quickly, therefore the 
terms with n = 1 mainly contribute to the reaction proba- 
bility. So, at frequencies where the probability gains maxima 
(f2' = 2m) contributions of those terms equal 



Or, 



,(0)| 



,sinh 2 (rY) 
T 2 r 2 



3j(h') 



Contributions of the terms with n = 2 



\CFr, 



2 sinh (Tr) + sin (6mu) g r) 2 

r 2 (r 2 + 36m 2 o;2) J2 ^ j 



obviously, are more than order of value smaller, in the case of 
uj g > T, i.e., when it makes sense to examine the interference 
in general. Thus, in order to make approximate assessments 
we omit the terms with n > 1. Then, for the same reason, for 
the ground state m = only contributions of the terms with 
n = are essential. It is those terms that make the contribu- 
tion independently of an AC MF: 



So = koo(0)| : 



; sinh 2 (rY) 



T 2 r 2 



As well, at a fixed frequency Q' = 2m* only terms with 
to = —to' = to* are essential in their contribution. Now the 
relative change of the concentration of the reaction products is 
easy to estimate at the MF frequency, e.g., f2' = 2m. Making 
note of J 2 _ 1 (/i') = J 2 (/i') and allowing for S = So + S m in 
this case, from (15) we arrive at 



,(0) 



-i -l 



(0) 



3l(h') 



As is seen, the magnitude of the magnetic effect depends on 
the ratio of the density matrix elements at the initial moment 
of time just after a gyro appears. For example, if the ground 
state and the state to (out of Zeeman's splitting) equipopulated 
at t = +0, then 



p = l- 



i + 3i(h<) 



(18) 



This function is shown in the Fig. 1, solid line. We conclude 
that the positions of the maxima of the amplitude spectrum of 
the magnetic effect do not depend (and the relative magnitude 
of the effect do) on the distribution of the initial populations 
of the gyro levels. 

The spectrum (17) determines only possible locations of ex- 
trema. A real form of the spectrum depends on the initial con- 
ditions for the density matrix, i.e., on the populations of levels 
of different rotational quantum number to. 

It is instructive to note that the molecule need not have 
a dipole moment J^i Qi r i f° r tne magnetic effect to appear. 
Rather, it is important that the 'charge moment of inertia' Q 



(10) be other than zero. This can be the case in the absence of 
dipole moment, e.g. for ionic rather than zwitterionic form of 
the molecule. 

The main properties of the gyro interference are identical 
with those of the ion interference, namely, (i) multiple peaks 
in the amplitude and frequency spectra, (ii) dependence of the 
positions of frequency peaks on the DC MF intensity, and (iii) 
independence of the positions of amplitude maxima on the AC 
MF frequency. 

We note that the interference of a molecular gyro has some 
features that differ it from the ion interference. Firstly, the 
peak frequencies are defined with respect to the gyral fre- 
quency uj g — a rotation equivalent of cyclotron frequency. 
Peak frequencies depend on the distribution of electric charges 
over the molecule and may deviate from harmonics and sub- 
harmonics of the cyclotron frequency. Secondly, the gyro 
rotation axis is fixed with respect to the shell, which intro- 
duces, in the general case, one more averaging parameter in 
the model. However, these features are not of principal signif- 
icance. The specific properties of the interference can always 
be calculated for any configuration of magnetic and electric 
fields, for rotation of biological systems and macromolecules 
involved, etc. 

There is the crucial feature of the gyro interference: molec- 
ular gyros are relatively immune to thermal shaking and may 
be effective biophysical targets for external MFs. 

As is seen from (16) the absolute magnitude of the mag- 
netic effect, where the latter is maximized by the MF param- 
eters, depends mainly on the value r) = Tt, which should be 
minimized for greater effects. The protein reaction time r and 
the MF frequency fl have to fulfill the relation VLt > 1 in or- 
der to manifest an interference. This and the properties of the 
function sinh 2 77/ 77 s lead to the condition of observability 



r- 1 > n- 



0.01s 



(19) 



for the ELF range. The following section examines if the con- 
dition is real. 



IV. ESTIMATING RELAXATION TIME FROM 
MOLECULAR DYNAMICS 

Computer simulation of molecular gyro behavior indicates 
that, for relaxation times of order 0.01 s, the size of cavity 
should be below 30 A. 

We consider the amino acid residue Phenilalanin, (Phe) 
C Q CgH5, as a gyro and look at the revolution of its benzene 
ring CeHs about the valence bond C a — Cp — see Fig. 4. This 
revolution may be thought of as a rotation in one plane of two 
rigidly bound point masses to = 26 m p (m p is the mass of 
proton), spaced by a = 2.42 A from one another, about their 
common center of gravity. 

We model the cavity by four heavy particles of mass M > 
m placed in the corners of a square (diagonal b > a) centered 
on the gyro axis, as shown in Fig. 5. We assume that these 
particles oscillate in the gyro rotation plane xy. Each particle 
moves in the potential well U(x i: j/j), where x t , yi is the de- 
viation of particle i from its equilibrium state. The Hamilton 
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FIG. 4: Schematic representation of a molecular gyroscope CaCeHs 
reduced to a model of a two-particle rotor of diameter a. 



FIG. 5: A two-dimensional model of a gyro in a molecular cavity of 
diameter b formed by four heavy particles of mass M. 



function for this system has the form 

4 



H = \lj> 2 + Y J [\M{i? i + yf) + V(4>,x i ,y i ) 

i=l 

+ U(xi,yi)}, 



(20) 



where / = ^ma 2 is the gyro moment of inertia, and cf> is its 
revolution angle. 

We take the potential of interaction of particle i with the 
gyro as the sum of two Lennard-Jones potentials 

VfaxtM) = e{[(r /n) 6 - l] 2 + [(r /r 2 ) 6 - l] 2 }, 

where ro is the equilibrium arm between a heavy and a light 
particle, n is the instantaneous distance of a heavy particle 
i to the first particle of the gyro, and r 2 is the distance to 
the second particle. The interaction of carbon atoms in poly- 
meric macromolecules is commonly described by Lennard- 
Jones potentials of the form 

VLj(r) = 4e [(a/r) 12 - (a/rf] 

with o- = 3.8 Aand e = 0.4937kJ/mol [27, 28]. Recognizing 
that each particle of the gyro consists of two carbon atoms, we 
let e = 1 kJ/mol ss 2e and r = 4.5 Aw 2 1 / 6 cr. 

The carrier potential for each heavy particle will be taken 
in the form 

tt/ ^ 1 x ' 2 + y 2 

U{x,y) = -K- , 

2 1 - (x l + y 2 )/Ro 

where K is the rigidity in particle-carrier interaction, and i? 
is the maximum possible deviation radius of a heavy particle. 



In a protein macromolecule, the rigidity of atomic displace- 
ments is K = 4 N/m. We consider two maximum displace- 
ment values: Rq = 1 A and Rq = oo. 

Assuming that heavy particles alone are connected with the 
thermostat, we obtain the equations of motion in the form 



Mjji 



9H 

dxi 
dyi 



■ T I Mx l + 

r r Myi + rji, 
i = 1,2,3,4, 



(21) 



where the system's Hamilton function is given by equation 
(20); £i and r\i are random normally distributed forces (white 
noise) describing the interaction of a heavy particle i with the 
thermostat, T T = 1/ t r is the friction factor, and t T is the parti- 
cle velocity relaxation time. The correlation functions of ran- 
dom forces are 



(m(ti)rij(t2)) 
(6(ti)%(ta)) 



2MT T k B T6 ij S(h - 1 2 ), 
2MT r k B TS ij 5(t 1 - 1 2 ), 
. 



Here, k B is the Boltzmann constant, and T is the thermostat 
temperature. 

We integrate the equation system (21) by the Runge-Kutta 
method to the fourth order of accuracy with a constant inte- 
gration step At. In this computation, the delta function S(t) 
is for \t\ > At/2 and 1/At for \t\ < At/2, that is, the in- 
tegration step corresponds to the correlation time of random 
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t (ps) 

FIG. 6: Current mean temperature of the molecular thermostat T\ 
(curve 1) and temperatures of its molecular neighborhood T2 (curve 
2) as functions of time. Thermostat temperature T = 300 K, cavity 
diameter b = 11 A, M = m, and Ro — 1 A. 

force. Therefore, to use a system of Langevin equations, we 
need that At <C t x . Let the relaxation time be t r = 0.2 ps, and 
the numerical integration step be At = 0.0025 ps. 

Let in the initial moment of time t = the system be in the 
fundamental state 

0(0) = 0o, Xi(0)=Ui, yi(0) = «i, (22) 
0(0) =0, ±i(0)=0, 2/i(0) = 0, * = 1,2,3,4, 

where the coordinates of a steady state, U > { u i, v i}f = i^ are 
determined as solutions to the minimization problem 

min : <f. >= 0, x\ = 0, 1/4 = . (23) 

4>,xi,...,yi 

Thus, at time zero, the molecular gyro is not thermalized. 

Our objective is to estimate the average time of gyro ther- 
malization. It corresponds to the relaxation time of gyro rota- 
tion in a thermalized molecular system. For this purpose, we 
numerically integrate the equations of motion (21) subject to 
the initial condition (22). 

The gyro thermalization at time t is characterized by its cur- 
rent temperature 

T 1 (t)=I{j> 2 (t))/k B , 

where brackets (■) imply averaging over independent realiza- 
tions of random forces £i(t), iji(t), i = 1,2,3,4. To obtain 
the average value, the system (21) was integrated more than 
10000 times. 

In turn, the thermalization of the system of heavy particles 
is characterized by its current temperature 

4 

B i=l 

The time dependence of these temperatures is presented in 
Fig. 6. At t = 0, the temperatures are Ti(0) = T 2 (0) = 0. 
Further on the time coordinate, they monotonously approach 
the thermostat temperature T = 300 K. 




10 15 20 25 30 35 

b (A) 

FIG. 7: Gyro thermalization time t± computed as a function of 
molecular cavity diameter b at M = m, Ro = 1 A (symbols 1); 
M — m, Ro — 00 (symbols 2); M = 100m, Ro = 00 (symbols 3), 
and extrapolation of this function to large b (curve 4). 



We will assume that the molecular subsystem is completely 
thermalized if its current temperature exceeds 0.99T. We 
determine the gyro thermalization time t\ as a solution of 
the equation T\{t) = 0.99T, and the time of heavy parti- 
cle system thermalization t2, as a solution of the equation 
Tiit) = 0.99T. The gyro is thermalized by interacting with 
the system of heavy particles, therefore its thermalization time 
will depend on the diameter b of the heavy particle system and 
will always exceed the heavy particle system thermalization 
time (t\ > £2)- Time t<i is almost independent of b and is 
dependent only on the relaxation time t T : t<z f=a 4i r . 

We analyzed the behavior of the system for Rq = 1 A, 00 
and M = m, 100m. The dependence of gyro thermalization 
time t\ on cavity diameter b is shown in Fig. 7. It is evident 
that, whatever the values of Ro and M, the thermalization 
time increases exponentially with b. If we extrapolate this de- 
pendence to the range of large b, we see that, at b = 28-32 A, 
the thermalization time, and hence the gyro relaxation time 
r _1 , will be of the order of seconds. With this size of cavity, 
the molecular gyro will revolve almost freely. 



V. CONCLUSION 

The molecular interfering gyroscope is a challenger for 
solving the kT-problem as a probable mechanism of magneto- 
biological effects. Indeed, the walls of a protein cavity do not 
interfere with the gyro degree of freedom directly via short- 
range chemical bonds. For cavities larger than 30 A in size, 
the contribution to the relaxation from the van der Waals elec- 
tromagnetic forces, induced by wall oscillations, is small. Ra- 
diation damping is negligibly small. Finally, the oscillations 
of gyro supports produce a zero moment of forces about the 
axis of rotation and do not affect the angular momentum. The 
gyro degree of freedom is very slow to thermalize, its dynamic 
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behavior is coherent, which gives rise to slow interference ef- 
fects. Of course, whether or not some more or less water-free 
cavities of the size of 30 A and larger do exist remains an open 
question, but, what is essential, ELF magnetic field bioeffects 
are no longer a paradox. 

The role of molecular gyros could probably be played by 
short sections of polypeptides and nucleic acids built inside 
globular proteins or in cavities between associated globules. 
In this respect it is interesting to look at Watson-Crick pairs of 
nitrous bases (adenine-thymine and guanine-cytosine) which 
bind DNA strands into a double helix as well as some other 
hydrogen-bound complexes of nitrous bases. Their rotations 
are hampered by steric factors. However, in the realm of activ- 
ity of special DNA enzymes, steric constraints may be lifted 
to allow a relatively free rotation of molecular complexes. It is 



not yet clear whether or not the gyro type of molecular struc- 
tures exists. They are unlikely to be detected by X-ray meth- 
ods since these require crystallization of proteins for structural 
analysis. In this state, the rotation would likely be frozen. 
Should a rotation be allowed, the mobile groups would not 
give clear cut reflections. Some other methods are needed that 
would work with native forms of proteins avoiding distortions 
due to crystallization. 

Generally speaking, the fact that the molecular gyro model 
gives a physically consistent explanation of MBEs proves in- 
directly its real grounds. Further studies should verify whether 
this conclusion is correct. In any case, today, the interfering 
molecular gyroscope is a single available mechanism to give 
explanations that would be physically transparent and gener- 
ally agreeable with experiments. 
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